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H^' Abstract 

^ . We address a number of unanswered questions in the M = (0,2) -deformed 

! CP(A^ - 1) model with twisted masses. In particular, we complete the program of 



X 
J3 



. solving CP(A^— 1) model with twisted masses in the large- limit. In hep-th/0512153 

nonsupersymmetric version of the model with the Zjv symmetric twisted masses was 
^ . analyzed in the framework of Witten's method. In larXivl0 ^803.0698 this analysis was 

O . extended: the large- iV solution of the heterotic A/" = (0, 2) CP(A^ — 1) model with 

no twisted masses was found. Here we solve this model with the twisted masses 
switched on. Dynamical scenarios at large and small m are studied (m is the twisted 
mass scale). We found three distinct phases and two phase transitions on the m 
plane. Two phases with the spontaneously broken Z^v-symmetry are separated by a 
phase with unbroken Z^. This latter phase is characterized by a unique vacuum and 
confinement of all U(l) charged fields ("quarks"). In the broken phases (one of them 
is at strong coupling) there are N degenerate vacua and no confinement, similarly 
to the situation in the A/" = (2, 2) model. Supersymmetry is spontaneously broken 
everywhere except a circle \m\ = A in the Zjv-unbroken phase. 

Related issues are considered. In particular, we discuss the mirror representation 
for the heterotic model in a certain limiting case. 
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1 Introduction 



Two-dimensional CP(A^ — 1) models with twisted masses emerged as effective low- 
energy theories on the worldsheet of non-Abelian strings in a class of four- dimensional 
M = 2 gauge theories with unequal (s) quark masses [H O |3l H] (for reviews see 
[5]). Deforming these models in various ways (i.e. breaking supersymmetry down 
to A/" = 1 and to nothing) one arrives at nonsupersymmetric or heterotic CP(A^ — 
1) modelsQ These two-dimensional models are very interesting on their own, since 
they exhibit nontrivial dynamics with or without phase transitions as one varies 
the twisted mass scale. In this paper we will present the large-A^ solution of the 
M = (0, 2) CP(iV — 1) model with twisted masses. As a warm up exercise we first 
analyze this model in the limit of vanishing heterotic deformation, i.e. the A^ = 
(2, 2) CP(A^ — 1) model with twisted masses (at A^ — )■ oo). The majority of results 
presented in this part of the paper are of a review nature and can be found in 
[2Zl EHl [ini I33]. In particular, Ref. [33] deals with the large- A^ limit in the Af = (2, 2) 
twisted mass deformed CP(A^ — 1) model. Our solution of the undeformed model 
exhibits two regimes - the strong coupling regime and the Higgs regime - with the 
crossover between them. We determine and briefly discuss the Argyres-Douglas 
points, an issue which was previously addressed in |33l [50]. We find it useful to 
collect the known results scattered in the literature, add some new nuances and, 
most of all, calibrate the basic tools to be exploited below, in the introductory part. 

Then we proceed to our main goal - the large-A^ solution of the heterotic de- 
formation. Both solutions (with and without deformation) that we present here are 
based on the method developed by Witten [6l [7] (see also [8]) and extended in [9] 
to include the heterotic deformation. For certain purposes we find it convenient to 
invoke the mirror representation (TUl [H] . 

An JV = (0,2)CP(A^ — 1) xC model on the string world sheet in the bulk 
theory deformed by fiA^ was suggested by Edalati and Tong [12|. It was de- 
rived from the bulk theory in [T3j (see also [TH [IS]). Finally, the heterotic JV = 
(0,2) CP(A^ — 1) model with twisted masses was formulated in [TB]. Its derivation 
from the microscopic bulk theory is under way [T7] . 

We report a number of exciting and quite unexpected results in the heterotically 

^ Strictly speaking, the full derivation of the heterotic CP(7V — 1) model with twisted masses, 
valid for arbitrary values of the deformation parameters, from the microscopic bulk theory, is still 
absent. 
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deformed CP(A^ — 1) model with twisted masses. The model has two adjustable 
parameters: one describing the strength of the heterotic deformation, and the other, 
m/A, sets the scale of the twisted masses. Dynamics of the model drastically changes 
as we vary the value of m, the parameter defining the twisted masses 



We discover three distinct phases on the m plane. In the first and the third phases, 
occurring at small and large values of |m|, the ZAr-symmetry of the model is sponta- 
neously broken. Correspondingly, there are N degenerate vacua and no confinement. 
In appearance, this is akin to what happens in the undeformed A/" = (2,2) model. 
However, the nature of these two phases is quite different. The first one occurs at 
strong coupling (small |m|) while the third one at weak coupling (large \m\). In 
essence, the latter is the Higgs phase. Surprisingly, at intermediate values of \m\ we 
find the Coulomb/confining phase, with unbroken ZAr-symmetry and unique vacuum. 
It is thoroughly investigated and the reasons for the photon to remain massless are re- 
vealed. Moreover, we find that at \m\ = A (in the the Coulomb/confining phase) the 
vacuum of the model is supersymmetric, while for all other values of m supersymme- 
try is spontaneously broken. At small and large values of the heterotic deformation 
parameter our solution is analytic. At intermediate values of this parameter it is 
semi-analytic: at certain stages we have to resort to numerical calculations. 

There are two phase transitions between the three distinct phases (Fig. [H]). We 
thoroughly analyze these phase transitions and argue that they are of the second 
kind. 

In addition to the large- A^ solution we address a number of related issues. For 
instance, in the limit of small heterotic deformation parameter we build the mirror 
representation for the J\f = (0, 2) CP(1) model. 

The organization of the paper is as follows. In a very short Section [2] we list 
general aspects of CP(A^ — 1) models. Section [3] introduces, in a brief form, our basic 
heterotic CP(A^ — 1) model with twisted masses in the gauged formulation most 
suitable for the large- A^ solution. (A discussion of the geometric formulation, which 
is also helpful in consideration of some aspects, is given in Appendix E). Section H] 
presents the large- A^ solution of the CP(A^ — 1) model with twisted masses and no 
heterotic deformation. Now, everything is ready for the comprehensive solution of 
the heterotic model. 




k = 0,1,2,... 
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In Sect. [5] we add a small heterotic deformation and analyze its impact on the 
large- solution. In Sect. [6] we find the analytic large- solution of the heterotic 
CP(A^ — 1) model with twisted masses in the limit of large deformations. Three 
phases and two phase transitions are identified. Intermediate values of the heterotic 
deformation parameter are studied semi-analytically and numerically. In Sect. [7] we 
focus on the second phase, namely, the Coulomb/confining regime. We explain here 
why, unlike two other phases, the photon remains massless. Section M is devoted to 
related issues and presents new results which are not necessarily based on large A^. 
Here we construct the mirror representation for the heterotic model with small values 
of the deformation parameter. Then we show, that unlike the undeformed model, the 
Veneziano-Yankielowicz-like and large- iV effective Lagrangians do not produce iden- 
tical results - a difference appears in the second order in the deformation parameter. 
Finally we present a new (albeit incomplete) derivation of the curve of marginal sta- 
bility in the large- limit (with Zn symmetric twisted masses). Appendices A and 
B explain the Euclidean vs. Minkowski notations and formulations. In Appendix 
C, for convenience of the reader, we compile and compare various definitions of the 
heterotic deformation parameters one can find in the literature. The large- A^ scaling 
laws for all these definitions are summarized here. In Appendix D we discuss global 
symmetries of the CP(A^ — 1) model with the Zjy-symmetric twisted masses. Finally, 
Appendix E gives a brief review of the geometric formulation of the heterotic model. 
Although it is only marginally used in this paper, it is convenient for related studies. 

Section M briefly summarizes our findings. 

2 Generalities 

A/" = (2, 2) supersymmetric CP(A^ — 1) sigma model was originally constructed [18] 
in terms of A/" = 1 superfields. Somewhat later it was realized [19] that = 
1 supersymmetry is automatically elevated up to A/" = 2 provided the target manifold 
of the sigma model in question is Kahlerian (for reviews see [201 EI]). The Witten 
index [22] of the CP(A^ — 1) model is A^, implying unbroken supersymmetry and A^ 
degenerate vacua at zero energy density. The CP(A^— 1) manifold is compact; there- 
fore, superpotential is impossible. One can introduce mass terms, however, through 
the twisted masses [23]. The model is asymptotically free [2l], a dynamical scale A 
is generated through dimensional transmutation. If the scale of the twisted masses 
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is much larger than A, the theory is at weak couphng. Otherwise it is at strong 
couphng. A priori, there are N distinct twisted mass parameters. However, in the 
absence of the heterotic deformation one of them is unobservable (see below). In this 
case the model is characterized by the coupling constant g^, the vacuum angle 6 and 
the twisted mass parameters mi, m2, ...,m7v with the constraint 

nil + 1712 + ••• + = . (2.1) 

By introducing a heterotic deformation, generally speaking, we eliminate the above 
constraint. The twisted masses are arbitrary complex parameters. Of special interest 
in some instances is the Z]sf symmetric choice 

(2'Kik\ 
, k = 0,l,2,...,N-l. (2.2) 

The set (12. 2p will be referred to as the ^Ar-symmetric masses. The model under con- 
sideration has axial 11(1)/? symmetry which is broken both by the chiral anomaly and 
the mass terms (see Apendix [D] for details). With the mass choice (12. 2p the discrete 
Z2N subgroup of this symmetry survives. We will see later that this symmetry is an 
important tool in studying phase transitions in the heterotic model. In analyzing 
some general aspects we will not limit ourselves to (12. 2p . The mass parameter m in 
Eq. (12. 2p can have an arbitrary phase. One can rotate away this phase at the price of 
generating a corresponding vacuum angle 6 (which effectively makes the dynamical 
scale parameter A complex. We will follow the convention in which A is kept real, 
while the phase is ascribed to m (in those issues where it is important). 

With the mass choice (12. 2 p the constraint (12. ip is automatically satisfied. Mostly 
in this paper we will consider the case of real and positive m. Sometimes however, 
we will relax this condition. 

Where necessary, we mark the bare coupling constant by the subscript and 
introduce the inverse parameter /3 as follows: 

/3o = \- (2.3) 

At large A^, in the 't Hooft limit, the parameter /3 scales as A^. 

There are two equivalent languages commonly used in the description of the 
CP(A^ — 1) model: the geometric language ascending to [19] (see also [21]), and the 
so-called gauged formulation ascending to [01 E]. Both have their convenient and 
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less convenient sides. We will discuss both formulations although construction of the 
expansion is more convenient within the framework of the gauged formulation. 
At |m|/A — 7- the elementary fields of the gauged formulation (they form an A^-plet) 
are in one-to-one correspondence with the kinks in the geometric formulation. The 
multiplicity of kinks - the fact they they enter in iV-plets - can be readily established 
[25] using the mirror representation [10]. We will discuss this in more detail later. 
We will review the gauged formulation of the model in the next section, while the 
geometric formulation is presented in Appendix [El 



3 The model 

In this section we will briefly review the gauged formulation of the model on which 
we will base the large- solution. An alternative geometric formulation useful for 
general purposes is presented in Appendix D. 



3.1 Gauged formulation, no heterotic deformation 

We start from the gauged formulation [6l [7] of the N = (2,2) CP(A^ — 1) model 
with twisted masses [23] setting the heterotic deformation coupling 7 = 0. This 
formulation is built on an A^-plet of complex scalar fields n* where i = 1,2,.. .,N. 
We impose the constraint 

nin' = 2(3. (3.1) 

This leaves us with 2A^ — 1 real bosonic degrees of freedom. To eliminate one extra 
degree of freedom we impose a local U(l) invariance n*(x) — ?■ e*°*-^^n*(x). To this end 
we introduce a gauge field which converts the partial derivative into the covariant 
one, 

d^^V,, = d^-iAf,. (3.2) 

The field is auxiliary; it enters in the Lagrangian without derivatives. The kinetic 
term of the n fields is 

C=\V^,n'\\ (3.3) 
The superpartner to the field ra* is an A^-plet of complex two-component spinor fields 

e, 

e = |^^ , (3.4) 
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subject to the constraint 

r^^e = 0, 6^' = 0. (3.5) 

Needless to say, the auxihary field has a complex scalar superpartner a and a 
two-component complex spinor superpartner A; both enter without derivatives. The 
full M = (2,2) -symmetric Lagrangian is0 

+ |V^nf + e,^a'^V^e + 2 5^ 



mi 

a — 

V2 



2 



(3.6) 



where are twisted mass parameters, and the limit Cg — )■ oo is implied. Moreover, 

a'^ = {l, zaa}, (3.7) 

see Appendix Rl 

It is clearly seen that the auxiliary field a enters in (13. 6 p only through the com- 
bination 

.--J.. (3.8) 

By an appropriate shift of a one can always redefine the twisted mass parameters in 
such a way that the constraint (12. ip is satisfied. The U(l) gauge symmetry is built 
in. This symmetry eliminates one bosonic degree of freedom, leaving us with IN — 2 
dynamical bosonic degrees of freedom inherent to CP(A^ — 1) model. 



3.2 Switching on the heterotic deformation 

The general formulation oiM = (0, 2) gauge theories in two dimensions was addressed 
by Witten in [7], see also [29]. In order to deform the CP(A^ — 1) model breaking 

■^This is, obviously, the Euclidean version. 
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M = (2, 2) down to A/" = (0, 2) we must introduce a right-handed spinor field 
(with a bosonic superpartner J^) whose target space is C, which is coupled to other 
fields as follows [El [13]: 

AC = U^dLCR + ^J" 

+ 2iuj\L(ji + 2iijj(jiXL — 2iuJ-'a — 2iuj^a, (3.9) 

where w is a deformation parameteilf]. 

This term must be added to the J\f = (2, 2) Lagrangian (13.61) . It is quite obvious 
that the dependence on (13.81) is gone. The deformation term (13. 9p has a separate 
dependence on a, not reducible to the combination (13. 8p . Therefore, for a generic 
choice, all N twisted mass parameters mi, m2, ...,mN become observable, Eq. (12. ip 
is no longer valid. 

Eliminating J-", and A, A we get 

AC = 4:\uj\'^\a\^ , (3.10) 
while the constraints (13. 5p are replaced by 

Ui^i = 0, ^^-n^ = 0, 

UiCR = V2ooCn, Un' = V2uU. (3.11) 

We still have to discuss dependence of the deformation parameter u. We want 
to single out appropriate powers of so that the large- iV limit will be smooth. From 
(I3.10p it is clear that u scales as \/N. 

One can restore the original form of the constraints (13. 5p by shifting the C,r fields, 
namely, 

^'r = U - j^V2ujnCR, = - l-y/2u:nCR- (3-12) 

This obviously changes the normalization of the kinetic term for C,ji, which we can 
bring back to its canonic form by a rescaling Cr., 

Ci?^(l-|7nCK, (3.13) 



■^Thc reader is referred to Appendix [C] where definitions of hetcrotic deformation parameters 
useful in different regimes are brought together. In this paper, only parameter w will be essentially 
exploited. 
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where the relation between 7 and uj as well as their A^- dependence are given in 
Appendix O As a result of these transformations, the following Lagrangian emerges 



+ 



+ ^ <y{ndknf - {nidRn)^L^L - {nidLn)^R^R 

+ (1 - l^fl^UL^R^R^L - ^R^R^L^L + iTpeieiCflC 



i7r) 



E 



m \n 



-im'\nf{^R^L)-tm'\nr{^LCR) 



I l„'|2/ 



(3.14) 



The sums over / above run from / = 1 to A^. If the masses are chosen Z^- 
symmetrically, see fl2.2p . this Lagrangian is explicitly Z^r-symmetric, see Appendix [Dl 
If all mi are zero, the model fl3.14p reduces to the Af = (0, 2) CP(A'" — 1) model 
derived in [13], see (13. 6p . Later on we will examine other special choices for the the 
mass terms. Here we will only note that with all mi ^ the masses of the boson 
and fermion excitations following from (I3.14p split [16]. Say, in the lo-th vacuum 



(0 



ferm 



M 



(0 

bos 



ferm 



l7l 



^,14 



1,2,...,A^; l^lo. (3.15) 

The splitting between the boson and fermion masses shows that A/" = (0,2) super- 
symmetry is spontaneously broken, see [TB] for further details. 

The model (I3.14p still contains redundant fields. In particular, there are A^ 
bosonic fields and A^ fermionic whereas the number of physical degrees of 
freedom is 2 x (A^ — 1)- One can readily eliminate the redundant fields, say, and 
, by exploiting the constraints (13. 5p . Then explicit Z^r-symmetry will be lost, of 
course. It will survive as an implicit symmetry. 
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4 Large- AT solution of the CP (AT — 1) model with 
twisted masses 



In this section we present the large- solution of the M = (2, 2) supersymmetric 
CP(A^ — 1) model with twisted masses fl3.6p . We consider a special case of mass 
deformation (12. 2p preserving the Zjv symmetry of the model. The = (2, 2) model 
with the vanishing twisted masses, as well as nonsupersymmetric CP(A^ — 1) model, 
were solved by Witten in the large- limit [6]. The same method was used in |30] 
to study nonsupersymmetric CP(A^ — 1) model with twisted mass. In this section 
we will generalize this analysis to solve the A/" = (2, 2) theory with twisted masses 
included. 

As was mentioned in the Introduction, many issues discussed in this section 
were previously addressed in [STJ [28l [TOl [33]. The large- limit in the Af = (2,2) 
supersymmetric CP(A^ — 1) model with twisted masses was treated in [33] • Moreover, 
the large- iV expansion, in fact, is not the only method in the studies of the A/" = (2, 2) 
supersymmetric CP(A^ — 1) model. Indeed, in this model exact superpotentials of the 
Veneziano-Yankielowicz type are known for arbitrary [351 IMl d EZl [28] . We use 
the large- A^ expansion in this section to prepare tools we will exploit later to solve 
the Af = (0, 2) supersymmetric CP(A^ — 1) model (for which no exact superpotentials 
are known). 

First let us very briefly review the physics of nonsupersymmetric CP(A^— 1) model 
revealed by the large- A^ solution [30]. In the limit of vanishing masses, the CP(A^ — 1) 
model is known to be a strongly coupled asymptotically free field theory [24]. A dy- 
namical scale A is generated as a result of dimensional transmutation. At large A^ it 
can be solved by virtue of the 1/A^ expansion [H]. The solution exhibits a "composite 
massless photon" coupled to A^ quanta n*, each with charge l/\fN with respect to 
this photon. In two dimensions the corresponding Coulomb potential is long-range. 
It causes linear confinement, so that only the nn pairs show up in the spectrum 
[Sn [B]. This is the reason why we will refer to this phase as "Coulomb/confining." 
In the Coulomb/confining phase the vacuum is unique and the Z^r symmetry is 
unbroken. 

On the other hand, if the mass deformation parameter m is ^ A, the model is 
at weak coupling, the field n develops a vacuum expectation value (VEV), there are 
A^ physically equivalent vacua, in each of which the Zn symmetry is spontaneously 
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broken. We refer to this regime as the Higgs phase. 

In Ref. [32] it was argued that (nonsupersymmetric) twisted mass deformed 
CP(A^ — 1) model undergoes a phase transition when the value of the mass parameter 
is ~ A, to the Higgs phase with the broken Zjy symmetry. In [30] this result was 
confirmed by the explicit large- solution. (Previously the issue of two phases and 
phase transitions in related models was addressed by Ferrari [331 El]-) 

In the M = (2,2) supersymmetric CP(A^ — 1) model, generally speaking, we do 
not expect a phase transition in the twisted mass to occur. In this section we confirm 
this expectation demonstrating that the Z^r symmetry is broken at all values of the 
twisted mass. (See, however, the end of Sect. 14. 2[ ) Still, the theory has two distinct 
regimes, the Higgs regime at large m and the strong-coupling one at small m0 

Since the action f l3.6p is quadratic in the fields and we can integrate over 
these fields and then minimize the resulting effective action with respect to the fields 
from the gauge multiplet. The large- limit ensures the corrections to the saddle 
point approximation to be small. In fact, this procedure boils down to calculating a 
small set of one-loop graphs with the rf and ^* fields propagating in loops. 

In the Higgs regime the field develops a VEV. One can always choose io = 
and denote = n. The field n, along with a, are our order parameters that 
distinguish between the strong coupling and Higgs regimes. These parameters show 
a rather dramatic crossover behavior when we move from one regime to another. 

Therefore, we do not want to integrate over n a priori. Instead, we will stick to 
the following strategy: we integrate over A^ — 1 fields n* with i ^ 0. The resulting 
effective action is to be considered as a functional of = n, D and a. To find the 
vacuum configuration, we will then minimize the effective action with respect to n, 
D and a. 

The fields n* and ^* (i = 1, ... A^ — 1) enter the Lagrangian quadratically, 



''At finite N there is no phase transition between these regimes. Instead, one has a crossover. 
This is explained after Eq. (|4.24|) . 




+ 




(4.1) 
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where the elhpses denote terms which contain neither n nor ^ fields. Hence, integra- 
tion over n* and C,^ in (13. 6p yields the following ratio of the determinants: 



nlT^et 



nii^det 



iD + \^/2a - 



rrii 



(4.2) 



where we dropped the gauge field Aj. which is irrelevant for the following determi- 



nation of vacuum structure 



The determinant in the denominator comes from the 



boson loops while that in the numerator from the fermion loops. Note, that the 



mass squared is given by iD + \V2 a 



while that of fermions is | \/2c" — 



VTii 



If supersymmetry is unbroken (i.e. = 0) these masses are equal, and the ratio of 
the determinants reduces to unity, as it should be, of course. 

Calculation of the determinants in Eq. (14. 2 p is straightforward. We easily get the 
following contribution to the effective Lagrangian: 



i=i I 



rrii 



In 



iD + I \/2a — rrii 



- \V2a- 



rrii 



In 



\V2a 



rrii 



(4.3) 



where quadratically divergent contributions from bosons and fermions do not depend 
on D and a and cancel each other. Here Muv is an ultraviolet (UV) cutoff. The bare 
coupling constant 2/3o in (13. 6p can be parametrized as 



N 



(4.4) 



Substituting this expression in (13. 6p and adding the one-loop correction (14. 3 p we 
see that the term proportional to iD In M^^ is canceled, and the effective action is 
expressed in terms of the renormalized coupling constant. 



1 ^\ iD+\V2a 
— > In 

4:71 ^ 



rrii 



A2 



(4.5) 



^Needless to say, this field is important in the spectrum calculation. 
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Assembling all contributions together and dropping the gaugino fields A we get 
the effective potential as a function oi n, D and a fields in the form 



Vcs = J ^iD + | -\/2cr — mo | ^ j |n 



2 



^ ^ iD + \ v^cr — rrii 



A2 



i=l 

1 \r^i r- 1 2 V^cr — mJ 1 

47r I A2 47r 



(4.6) 



Now, to find the vacua, we must minimize the effective potential (14. 6p with respect 
to n, D and a. In this way we arrive at the set of the vacuum equations, 

\n\' = 2/3rcn, (4.7) 

\^iD + \V2a - mo\^^ n = 0, (4.8) 

1 r- ^D + \V2a '''^ 



(v^a-mo)|nr--5^(v^a-m^)ln^L I = q 



(4.9) 



where 2/3ren is determined by Eq. (14.51) . 

From Eq. (14. 8 P it is obvious that there are two options: either 



D + \V2a-mo\ = (4.10) 



or 

n = 0. (4.11) 

These two distinct solutions correspond to the Higgs and the strong-coupling regimes 
of the theory, respectively. Equations (I4.7p - (l4.9p represent our master set which 
determines the vacua of the theory. 
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4.1 The Higgs regime 

Consider first the Higgs regime. For large m we have the solution 

D = 0, V2a = mo, \n\^ = 2f3,^^. (4.12) 

The first condition here, D = 0, means that A/" = (2, 2) supersymmetry is not broken 
and the vacuum energy is zero. Integrating over n's and ^'s we fixed rP = n. Clearly, 
alternatively we could have fixed any other Then, instead of (I4.12p . we would 

get 

= 0, ^/2a = m,„, In'op = 2/3,en , (4.13) 

demonstrating the presence of degenerate vacua. The discrete chiral Z2N symme- 
try (1D.6P is broken by these VEV's down to Z2. Substituting the above expressions 
for D and a in (14. 5 P we get the renormalized coupling 

1 I 12 AT 

1 \ ^ , rrin — mA i\ ^ m 

i=l 

where we calculated the sum over i in the large- limit for the special choice of 
masses (12. 2p . 

In each vacuum there are 2(A^ — 1) elementary excitationsl^ with the physical 
masses 

Mi = |mj-mij, iy^io- (4.15) 

In addition to the elementary excitations, there are kinks (domain "walls" which are 
particles in two dimensions) interpolating between these vacua. Their masses scale 
as 

The kinks are much heavier than elementary excitations at weak coupling. Note 
that they have nothing to do with Witten's n solitons [6] identified as the fields 
at strong coupling, see Sect. 18.31 

Since jn*"!^ = 2/3rcn is positively defined we see that the crossover point is at 
m = A. Below this point, the VEV of the n field vanishes, and we are in the strong 
coupling regime. 



^Hcrc wc count real degrees of freedom. The action p.6p contains N complex fields n'. The 
phase of can be eliminated from the very beginning. The condition = 2/3 eliminates one 
extra field. 
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4.2 The strong coupling regime 

For small m the solutions of Eqs. ( \A.7\} - (14. 9p can be readily found, 

D = 0, n = , 2/3,en = ^ ^ a5 ^ = 0- (4-17) 

i=l 

Much in the same way as in the Higgs regime, the condition D = means that 
A/" = (2, 2) supersymmetry remains unbroken. 

Note that at large A^, the summation in f l4.17p can be extended to include the 
i = term, 

2/5- = 4^ E 1^ = ' (4-18) 

because (as we will show below) \/2a ~ A in this regime and is not close to any of 
rrii at \m\ < A. 

The last equation can be identically rewritten as 

Yl |v^fr-m, = A^. (4.19) 

1=0 

For the Zjy-symmetric masses Eq. (14.191) can be solved. Say, for even A^ one can 
rewrite this equation in the form 



V2a 



N 



A^ , (4.20) 



due to the fact that with the masses given in (12. 2p 

^ mi = , 

mirrij = , 

mi^mi^...mi^_^ = 0, (ii ^ i2 ^ ■■■ ^ iN-i) ■ (4.21) 

«l,22vi«iV-l 

Equation (14.200 has A^ solutions 

v^a = (A^ + m^)'/^exp(^^^, A; = 0, A^ - 1, (4.22) 
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n 




A 



Figure 1: Plots of n and a VEVs (thick lines) vs. m in the TV = (2, 2) CP(A^ - 1) model 
with twisted masses as in (j2.2p . 

where we assumed for simplicity that m = rriQ is real and positive. (This is by no 
means necessary; we will relax this assumption at the end of this section.) Note that 
the phase factor of a in (14.221) does not follow from (I4.19p . Rather, its emergence 
is explained by explicit breaking of the axial U{1)r symmetry down to Z2N through 
the anomaly and non-zero masses (12.21) . see Appendix [Dl with the subsequent spon- 
taneous breaking of Z2N down to Z2. Once we have one solution to ( 14.19^ with the 
nonvanishing a we can generate all solutions (I4.22p by the Z2N transformation [6] . 

Although we derived Eq. (I4.19P in the large- A^ approximation, the complexified 
version of this equation. 



is in fact, exact, since this equation as well as the solution (I4.22p follow from the 
Veneziano-Yankielowicz-type effective Lagrangian exactly derived in the A/" = (2, 2) 
CP(A^ — 1) model in |35l [361 El EZl 128]. The Veneziano-Yankielowicz Lagrangian 
implies (I4.23P even at finite A^. 




(4.23) 



i=0 
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The solution (14.22^ shows the presence of N degenerate vacua. Since a 
in all these vacua the discrete chiral Z2N symmetry is broken down to Z2 in the 
strong-coupling regime, much in the same way as in the Higgs regime. This should 
be contrasted with the large- iV solution of the nonsupersymmetric massive CP(A^ — 
1) model [30]. In the latter case, a = in the strong coupling phase, therefore, the 
theory has a single vacuum state in which the Z2N symmetry is restored. This is a 
signal of a phase transition separating the Higgs and Coulomb/confining phases in 
the nonsupersymmetric massive CP(A^ — 1) model [30] . 

In fact, in the large- approximation the formula (I4.22p can be rewritten as 

r- f2-Kik\ f A, m < A , , 

V2a = exp — - X A; = 0, AT _ 1 (4.24) 

\ JM J [ m, m> A 

with the exponential accuracy O (e^^) . Note that at large m this formula reproduces 
our result (14.131) obtained in the Higgs regime. In the limit m — )■ it gives Witten's 
result [6]. 

The VEVs n and a as functions of m are plotted in Fig. [1] These plots suggest 
that we have discontinuities in derivatives over m for both order parameters. Taken at 
its face value, this would signal a phase transition, of course. We note, however, that 
the exact formula (I4.22p shows a smooth behavior in a. Therefore, we interpret the 
discontinuity in (14.241) as an artifact of the large- A^ approximation. The crossover 
transition between the two regimes becomes exceedingly more pronounced as we 
increase A^ and turns into the second-order phase transition in the limit N ^ 00. 
We stress again that the Z2N symmetry is broken down to Z2 in the both regimes. 

There is one interesting special point in Eq. (14.231) . Relaxing the requirement of 
reality of the parameter mo we can choose the product 

n(-m,) = A^. (4.25) 

i=0 

At this particular point Eq. (I4.23P reduces to cr^ = with the solution 

a = 0. (4.26) 

All A^ vacua coalesce! This is a two-dimensional "reflection" of the four-dimensional 
Argyres-Douglas (AD) point [38| [39]. 
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4.3 Generic twisted masses and the Argyres— Douglas points 

In this subsection we briefly describe the AD points in the undeformed M = (2,2) 
CP(iV — 1) model. At these points one or more kinks interpolating between different 
vacua of the model become massless. These points determine a nontrivial conformal 
regime in the theory. We use the complexified version of the vacuum equation (14.23^ 
appropriate for a generic choice of the twisted mass parameters. 

Let us have a closer look at this equation given a set of arbitrary masses. Our 
task is to find the values of mass parameters such that two roots of this equation 
coalesce, (Xi = cr2- Near the common value of a Eq. fl4.23p can be simplified, namely, 

{V2a-m,,) = ^ ^ (4.27) 

where ^ 

"ii2 = -{mi+m2), Ami2 = nii - m2 (4.28) 

Equation f l4.27p gives 



v^ai,2 = mi2 ± + Ae\. (4.29) 



4 



Two vacua coalesce if the square root vanishes. 



- Am?2 JJ (mi2 - mi) = 4A^ . (4.30) 

At this AD point one of kinks interpolating between the vacua at ai and a2 
becomes massless. 

Similarly, one can consider more complicated AD points in which more than two 
vacua coalesce. At these AD points more than one kink becomes massless. The 
point f l4.25p corresponds to a very special regime in which all vacua coalesce (for 
the ZTv-symmetric choice of masses on the circle (12.20 ). At this point in the mass 
parameter space one of kinks interpolating between each two "neighboring" vacua 
becomes massless. This AD point was studied previously in [SD] . We remind that the 
A/" = (2, 2) supersymmetric CP(A^ — 1) model is an effective theory on the world sheet 
of the non-Abelian string in A/" = 2 SQCD (with the U(A^) gauge group and Nf = N 
flavors P, 121 El H] ) . Therefore, the massless kinks at the AD points in two dimensions 
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correspond to massless confined monopoles at the AD points in four- dimensional bulk 
theory. 

We pause here to make a remark unrelated to the Argyres-Douglas points. As- 
sume one has a (nearly) generic set of twisted masses subject to a single constraint 

N-l 

n(-m.) = A^. (4.31) 

Then Eq. fl4.23p has a solution cr = 0, with other N — 1 solutions a ^ 0. Now, if we in- 
troduce the heterotic deformation ~ a^, the vacuum a = remains supersymmetric, 
while in all other vacua supersymmetry is broken. 



5 CP(A^ — 1) model at small heterotic deformations 

Now, we switch on the heterotic deformation which breaks A/" = (2, 2) supersym- 
metry down to A/" = (0,2). In this section we will assume this deformation to be 
small limiting ourselves to the lowest nontrivial order in the heterotic deformation. 
All preparatory work was carried out in Sect. |H Therefore, here we can focus on the 
impact of the heterotic deformation per se. 

To determine the effective action allowing us to explore the vacuum structure of 
the heterotic model, just as in Sect. HI we integrate over all but one given n' field 
(and its superpartner ^'). One can always choose this fixed (unintegrated) field to 
be = n. Assuming a and D to be constant background fields, and evaluating the 
determinants one arrives at the following effective potential (see Eq. fIS.lOp ): 

Ves = / (fxl (iD + |V2cr-mon 



1 ^— ! / , r- •i2\ iD + \\f2o — mM^ , ^ 

- 5^ (.D + I V2a - mf ) In L (5.I) 

i=l 



1 r- -12 \\f2o-m^\ \ , , N , 

+ — y\V2a- mf In ^- L + w {N - 1) + —-ula 

47r I A2 47r 27r ' 

i=l 



where we have introduced a deformation parameter 

Sir 
iV 



Sir , ,^ / N 

u = — (5.2) 
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Note that although \u\'^ grows as 0{N) for large A^, the parameter u does not scale 
with and so is more appropriate for the role of an expansion parameter. 

The above expression for Vcs replicates Eq. (14 .Gp except for the last term repre- 
senting the heterotic deformation. Now, to find the vacua, we must minimize the 
effective potential (15. ip with respect to n, D and a. The set of the vacuum equations 
is 



n 



^-SE'" 4^ L = o, (5.3) 



%D ^\\P1g -mX]n = ^ , (5.4) 



(V2a - mo) > (V2a - m') In ^ ^ + ^-u V2a = . 

(5.5) 

It is identical to the master set of Sect. H] with the exception of the last term in 
Eq. (15. 5p . Equation (15. 4p is the same; hence we have the same two options: either 

iD + \V2a - mo\^ = (5.6) 



or 



n = 0. (5.7) 



Since the deformation parameter is assumed to be small, we will solve these equations 
perturbatively, expanding in powers of u, 

n = n(°^ + u ■ n^^^ + . . . , 

iD = iD^^^ + u ■ iD^^^ + . . . , 

a = + u ■ a^^^ + ... . (5.8) 

Here n^°\ D^^) and a^") constitute the solution of the A/" = (2,2) CP(A^ - 1) sigma 
model, in particular D^^^ = in both cases (15. 6p and (15. 7p corresponding to the 
Higgs and the strong-coupling regimes of the theory, respectively. We remind that 
the mass parameters are chosen according to (12. 2p . 
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5.1 The Higgs regime 



The large- supersymmetric solution of the M = (2, 2) CP(A^ — 1) sigma model in 
the Higgs phase is given in Eqs. fl4.13p and f l4.14p . Expanding Eqs. (15. 3p - (15. 5p to 
the first order in u, we calculate 







\n 



(0)12 



n 



(1) 



2 m 



N 



77,(0) 1 77,(0) 1 2 327r2 ^-^ m — 



N-1 



With masses from (12. 2p we then obtain 



A^- 1 



i=l 



m — m'- 



2m 



N 
2m 



+ 0(11 



Using this, we simplify the solution (15. 9p . 



V2i 



a 



iD 



m 1 



m 



u/2 ' 
In m/ A ^ 



+ 



u/2 
Inm/A 



n 



2/3ren 

This is in the Higgs phase, where 

2/3(°) 

t^ren 



{\nm/Ay 



2n \AJ 



(5.9) 



(5.10) 



(5.11) 



Eqs. (15. lip agree with numerical calculations of the solution of the vacuum equations 
in the Higgs phase. 



5.2 Strong coupling 

Our starting point is the zeroth order in u solution (Sect. |4]), 

,(0) = 0, = 0, V2a(°) = A-e'"-^ 



(5.12) 
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where 



A 



k\i + O (e~^)] at 



oo . 



(5.13) 



At strong coupling n vanishes exactly, not only in the zeroth order in u. Omitting 
the details, the first order solution to the vacuum equations fl5.3p . fl5.5p is given by 
(in conjunction with n = 0) 



2^(0) 



N ^ 







N-l 



A^ ^ x/2a(o) - 
1=1 * 



AT Z^i=l 

+ h.c. 



(5.14) 



V2<f(0)- 



|V2<T(0)-m*|2 



We use the following relations to simplify the above formulas when masses are set as 
in (O: 



N-l 

-T — - 

k=o ^ " 



27rik 

— ae N 



a 



N 



N-l 

-T- 



1 



1 + a 



N 



1 



(5.15) 



A^ f-^ (l + a2) -2acos^ 1 - 1 - I - a' 

This enables us to present the results for m ^ A in the following quite simple form: 

n = 0, iD = uA"^ + ... , 

u A2 



2iTil 

Ae ^ 



1 



2 A2 - 777,2 



+ 



(5.16) 



We complement these formulas for the strong coupling phase by finding the approx- 
imate solutions now as expansions in tti^ parameter, assuming ttj to be small. We 
obtain 



v/2ct = e^A e-"/2 
iD = A2(1 - e"") 



777 

-^shM/2 



+ ... 



O 



777^ \ 

A4 J 



(5.17) 



where u does not need to be (too) small anymore. 
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Just a brief look at the Higgs phase solution ( IS.lip and the strong coupling 
phase solutions f l5.16p and fl5.17p reveals, that these expansions blow up when one 
approaches m ^ A! While the exact solutions are expected to be finite for all m, our 
approximations cannot be trusted at m = A. This is the first sign that something is 
going on at these values of masses. As we will later see from the large-w solutions, as 
well as from the numerical solution of the vacuum equations, the theory experiences 
a double phase transition as m goes from the area below A towards m 3> A. 

6 Heterotic CP(A^— 1) model at large deformations 

Now it is time to study equations fl5.3p - flS.Sp in the opposite limit of large values 
of the deformation parameter u ^ 1. We will see that our theory has three distinct 
phases separated by two phase transitions: 

(i) Strong coupling phase with the broken Z^r symmetry at small m; 

(ii) Coulomb/confining Z^r-symmetric phase at intermediate m (the coupling con- 
stant is strong in this phase as in the case (i)); 

(iii) Higgs phase at large m where the Zj^f symmetry is again broken. 

As previously, we assume that mass parameters are chosen in accordance with 

6.1 Strong coupling phase with broken Zn 

This phase occurs at very small masses, namely, 

m<Ae""/^ M>1. (6.1) 

In this phase we have 

\n\ = 0, 2/3,en = 4^ E 1^ ^^A^ = ^ ■ (6-2) 

1=1 

As we will see momentarily, a is exponentially small in this phase. Masses are also 
small. Then the second equation in (16. 2 p gives 

iD ^ A^ (6.3) 
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m''' = m e'' 




Figure 2: Electrostatic analog problem. The circle of radius m (see Eq. (j2.2p ) is homoge- 
neously populated by "electric charges," namely, Ok = 2-Kk/N where /c = 0, 1,2, — 1. 
We then must calculate the electrostatic potential at the point x. 

With this value of iD we can rewrite Eq. (15. 5p in the form 



N-l 



^ (y2a - In ^ = N (y2a^ u . (6.4) 

i=i I "\/2cr — mi I 

The following trick is very convenient for solving this equation. 

Let us consider an auxiliary problem from static electrodynamics in two dimen- 
sions. Assume we have N equal "electric charges" evenly distributed over the circle 
depicted in Fig. [21 In the limit of large one can consider this distribution to be 
continuous (and homogeneous). The task is to find the electrostatic potential at the 
point X on the plane. 

It is not difficult to calculate the potential of a charged circle of radius m centered 
at the origin in two-dimensional electrostatics. Representing x by a complex number 
we get (in the large- iV limit) 

1 ^i;^ I In \x\ > m 

-}_^ln |x-m,|2 = < . (6.5) 

Now, to obtain the left-hand side of (16. 4p we must integrate ( 16. 5p over x and then 
substitute x = \/2a. In this way we arrive at 



2^ i^V2a - m,j In 



|V2a-m,|2 |y2a(ln^-l), \V2a\ < m ' 

(6.6) 
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Outside the circle the potential is the same as that of the unit charge at the origin. 
Inside the circle the potential is constant. 

Substituting Eq. (16. 6 p in (16. 4 p at m < |a/2o"| (i.e. outside the circle) we get 

y2(^) = e^'- Ae-"/2, A; = 0, (A^ - 1). (6.7) 

The vacuum value of a is exponentially small at large u. The bound m < \V^(7\ 
translates into the condition (16. ip for m. 

We see that we have A^ degenerate vacua in this phase. The chiral Z2N symmetry 
is broken down to Z2, the order parameter is (a). Moreover, the absolute value of a 
in these vacua does not depend on m. In fact, this solution coincides with the one 
obtained in [9] for m = 0. This phase is quite similar to the strong coupling phase 
of the J\f = (2, 2) CP(A^ - 1) model, see fO^ . The difference is that the absolute 
value of a depends now on u and becomes exponentially small in the limit u ^ 1. 

The vacuum energy is positive (see Eq. (16. 3p ) - supersymmetry is broken. We 
will present a plot of the vacuum energy as a function of m below, in Sect. 16. 2[ 

6.2 Coulomb/confining phase 

Now we increase m above the bound (16. ip . From (16. 6p we see that the exponentially 
small solution to Eq. (16. 4 p no longer exist. The only solution is 

(a) = 0. (6.8) 

In addition, Eq. (16. 2 p implies 

\n\=0, iD = k^-m^. (6.9) 

This solution describes a single Z^r symmetric vacuum. All other vacua are lifted 
and become quasivacua (metastable at large A^). This phase is quite similar to the 
Coulomb/confining phase of nonsupersymmetric CP(A^ — 1) model without twisted 
masses [6] . The presence of small splittings between quasivacua produces a linear ris- 
ing confining potential between kinks that interpolate between, say, the true vacuum 
and the lowest quasivacuum [32], see also the review |1D]. As was already men- 
tioned, this linear potential was interpreted, long time ago [311 E] , as the Coulomb 
interaction, see the next section for a more detailed discussion. 

As soon as we have a phase with the broken symmetry at small m, and the 
^TV-symmetric phase at intermediate m the theory experiences a phase transition 
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A2 




m 



Figure 3: Vacuum energy density vs. m. The dashed line shows the behavior of the energy 
density ()6.10p extrapolated into the strong coupling region. 

that separates these phases. As a rule, one does not have phase transitions in super- 
symmetric theories. However, in the model at hand supersymmetry is badly broken 
(in fact, it is broken already at the classical level [H]); therefore, the emergence of a 
phase transition is not too surprising. 

We can calculate the vacuum energy explicitly to see the degree of supersymmetry 
breaking. Substituting (16. 8 p and (16. 9p in the effective potential (15. ip we get 



The behavior of the vacuum energy density -Evac vs. m is shown in Fig. |3l 

Evac is positive at generic values of m, as it should be in the case of the sponta- 
neous breaking of supersymmetry. Observe, however, that the vacuum energy density 
vanishes at m = A. This is a signal of A/" = (0, 2) supersymmetry restoration. To 
check that this is indeed the case - supersymmetry is dynamically restored at m = A 
- we can compare the masses of the bosons and their fermion superpartners 
From (13. 6 p we see that the difference of their masses reduces to iD. Now, Eq. (16. 9 p 
shows that iD vanishes exactly at m = A. 

This is a remarkable phenomenon: while A/" = (0,2) supersymmetry is broken at 
the classical level at m = A, it gets restored at the quantum level at this particular 



7?' 



vac 



Coulomb 




(6.10) 
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point in the parameter space. This observation is imphcit in |4T] where a Veneziano- 
Yankielowicz-type (VY-type) superpotential [12] for M = (2,2) CP(A^ — 1) model 
(see [351 ESI E]) was extrapolated to the Af = (0, 2) case. 

We pause here to make an explanatory remark regarding Fig. [3] and Eq. fl6.10p . 
The plot of E^^^^""^^ is presented in this fi gure assuming the parameter m to be real 
(we follow this assumption in the bulk of the paper). In fact, m can be viewed as a 
complex parameter, the phase of m being interpreted as a ^ angle. A straightforward 
examination shows that for the complex values of m Eq. f l6.10p must be replaced by 



r'Coulomb ^ 

An 



A^ — \m\^ + Imp In 



m 



2 



A2 



This means that the vacuum is supersymmetric (i.e. = 0) on the curve 

|mp = A^. 

Now we turn our attention to what happens with a at the strong/Coulomb phase 
transition point. More detail on that will be given in Section [6.41 with the help of 
numerical calculations, however, at large u the behavior of a can be analyzed just 
by inspecting Eq. (16. 4p . To solve this equation we can evaluate the sum in it using 
Eq. (16. 6p . In place of the massive quantities a, A and m it is convenient to introduce 
dimensionless variables 



S = = m/A. (6.11) 



A 



Then Eq. (16. 4p turns into 



+ In^^j . (6.12) 



Instead of solving for S one can use (I6.12p as a solution for m with respect to 
cr. Figure H] illustrates the dependence (16.120 . We now treat this graph as the 
dependence of a on m. In particular, for zero masses, \/2a = Ae~"/^, as it should 
be. Following from right to left, as the mass grows, a decreases, and at some point 
the curve bends downward. This happens at jj? = S'^ = e"*-^^'"'', as can be seen from 
Eq. (I6.12p . At this point the derivative da /dm becomes infinite, which is indicative 
of a phase transition. The behavior of a near this point is circle-like, 

V2a ~ Ae"^ + -Ay^e-("+i) - m? / . (6.13) 
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Figure 4: Dependence of versus a^. The rescaled variables as in text are = m/A 
and S = V^a/A. 

Equations fl6.12p and fl6.13p together with Fig. H] are approximate, but quahtatively 
they very closely demonstrate what happens to a in reality. Namely, a monotonically 
decreases with increase of mass, until it experiences a vertical bend, at which point 
it drops down to zero, the fact that cannot be seen from Eq. fl6.12p . Rather, this is 
seen in the Coulomb phase, and reproduced with numerical solution in Section 16. 4[ 
Despite the fact that apparently a experiences a drop, the energy density does not 
(see Fig. [6]), and therefore the phase transition is of the second order. 

6.3 Higgs phase 

The Higgs phase occurs in the model under consideration at large m. Below we will 
show that the model is in the Higgs phase at 

m > , if n > 1 . (6.14) 

In this phase \n\ develops a VEV. From Eq. (15. 4p we see that 

iD = - |v^cr-mo|^ (6.15) 

To begin with, let us examine Eqs. ( 15. 3p - (15. 5 p far to the right from the boundary 
(I6.14p . i.e. at m ^ ^/uA. In this regime we can drop the second logarithmic term 
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in (I5.5p . This will be confirmed shortly. The first term is much larger because it is 
proportional to /3ren which is large in the quasiclassical region (see Eqs. (14. 7p and 
fHl2ll ). Then Eq. (ES]) reduces to 



^ r 
+ — nV2o- 

47r 



0, 



implying, at large u 



V2i 



a 



Stt 



-.J /3rcn 1 



mo 



(6.16) 



(6.17) 



where we take into account that \a\ ^ m, the fact justified a posteriori. Equation 
fl6.17p applies to the k = vacuum. It is obvious that the solution for other — 1 
vacua can be obtained from fl6.17p by replacing mo — rrii^ where i^ = 1, {N — 1). 

Thus, we have degenerate vacua again. In each of them \a\ is small (~ m/u) 
but nonvanishing. The Z2N chiral symmetry is again broken down to Z2. Clearly, the 
Higgs phase is separated form the Coulomb / confining phase (where Z2N is unbroken) 
by a phase transition. 

To get the vacuum expectation value of we must analyze the logarithms in 
Eq. f l5.3p and f l5.5p with a better accuracy: a in the numerators cannot be neglected. 
We must keep the terms linear in a. Since the solution for a is real, see fl6.17p . we 
can rewrite the logarithm in (15. 3p as follows: 



\n(iD+ \V2i 



a 



m 



In 



4^2, 
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N 



k = 1,...,N -1. (6.18) 



where is the phase of m^, see Fig. [21 On the other hand, Eq. fl4.14p can be 
presented in the form 



1 ^""^ a 
— > In 4 sin^ ( — 



Thus, we conclude that 



k=l 



\n\ 
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(6.19) 
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N m^ 
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(6.20) 
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in each of the vacua in the Higgs phase. Here the last (rather rough) estimate 
follows from f l6.17p . 

Our next task is to get an equation for /3ren (en route, we will relax the constraint 
m ^ v^^)- To this end we must examine Eq. (15. 5p . including the logarithm into 
consideration. We will expand the numerator neglecting O(cr^) terms, while in the 
denominator we can set a = right away. Then the summation in (15. 5p can be 
readily performed using the formula 

— (mo — rrifc) In 4 sin^ ("~~) ~ ^ ~^ ' (^-21) 

fc=0 

which follows, in turn, from Eq. (16.61) . As a result, we arrive at 

V2au = m (3ron + 1^ • (6.22) 

The only approximation here is u ^ 1, plus, of course, Eq. (I6.14p . Combining 
Eqs. (I6.22p and (I6.20p we obtain the following relation for /3ren: 

|/3.en-ln(|/3.en + l)=ln^. (6.23) 

Strictly speaking, Eq. (I6.23P has two solutions at large m, deep inside the Higgs 
domain. The smaller solution corresponds to negative /3ren- Since |?7,p = 2/3ren is 
positively defined we keep only the larger one. At m ^ y/uA /3rcn is large and is 
given (with the logarithmic accuracy) by the last estimate in Eq. (I6.20p . As we reduce 
m, at m = y/uA, two solutions of (I6.23P coalesce. At smaller m they become complex. 
Thus m = y/uA is indeed the phase transition point to the Coulomb/confining phase. 
At this point /3rcn = 0, which coincides with its value in the Coulomb/confining 
phase, see (16. 9p . Thus, /3rcn is continuous at the point of the phase transition, while 
its derivative over m is discontinuous. 

Calculating the vacuum energy in this phase we get 



pHiggs _ jj_ 

in 



m In — + A^ - + O — In — 



A2 ' A^ ) ] 



(6.24) 



The vacuum energy in the Higgs phase is non-zero so the M = (0,2) supersymmetry 
is broken. It was observed earlier in [TB] on the classical level, see also fl3.15p . The 
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Figure 5: The dependence of a on m. The sohd curve shows the numerical solution of 
the equations ()5.3p and ()5.5p in the strong coupling phase (|n| = 0) for u = 2. The upper 
dashed curve shows the small-m limit ()5.17|) . The lower dashed curve refers to the large-u 
solution obtained as a numerical evaluation of Eq. ()6.12p . The dash-dot ^/2a = m line is 
given for reference. To the right of the drop, a is identically zero. 

vacuum energy density in all phases is displayed in Fig. El We can see that the 
expression for energy fl6.24p at large u coincides with Eq. fl6.10p . which signifies that 
the phase transition is of the second order. In Section 16.4.21 we will analyze the 
transition in more detail. 

6.4 Evaluation at arbitrary magnitude of deformation 

In this subsection we will grasp a picture of what happens to the order parameters 
in phase transition regions when the deformation u is not necessarily large or small. 
Although we will be able to acquire an insight into the solution of the vacuum equa- 
tions in the Higgs phase analytically, we will need to revoke numerical calculations to 
examine the energy density and the coupling constant. In the Strong phase, even less 
can be done analytically, and we will resort to numerical methods almost entirely. 
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Figure 6: Energy density versus m in the strong coupling and Coulomb phase for u = 2. 
The dashed curve displays the energy density of the Coulomb/confining phase, Eq. (I6.10p . 

6.4.1 Strong/Coulomb phase transition 

To see what happens at the strong-to-Coulomb phase transition we solve the vacuum 
equations (15. 3p - (15.51) numerically, with n = 0. Figure [5] shows the dependence of 
cr on m in the strong coupling phase. The main graph is compared to the small-m 
solution (15.171) and the large-w dependence (I6.12p (still, solved numerically). One 
convinces that the large-w solution given by Eqs. (16. 4p and (I6.12p describes the true 
solution really well. The curve for \/2a monotonically decreases from the value 
Ae~"/^, until it meets the line \/2a = m, at which point a turns down 90° and 
drops vertically to zero. This is the strong-Coulomb phase transition point, which is 
approximately located at m = \/2a ~ e'^^^"-*/^. To the right of the phase transition 
point a is identically zero. 

Figure [6] shows the energy density in the strong coupling and Coulomb/confining 
phase. The curve clearly displays that the energy does not experience any jumps at 
the phase transition point m ^ e-(i+")/2_ That is, the phase transition is at most of 
the second order. The curve, however, does apparently experience a break of incline 
at that point. To the right of the phase transition, the numerical curve exactly 
overlays the Coulomb phase energy density Eq. (I6.10p . 

We further clarify the position of the phase boundary numerically, by solving the 
vacuum equations with the condition m = \/2a. The phase boundary happens to 
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coincide exactly with the left branch of the solution of the equation 

u = i4 - In/i^ - 1, 

which will arise in the next subsection when we will be analyzing the Higgs phase 
transition. At that point the full phase picture of the theory will be explicated. 

6.4.2 Higgs phase 

We revoke the dimensionless variables S and fi introduced in Eq. ( 16. lip . The system 
of vacuum equations fl5.3p - fl5.5p . written in terms of these variables, becomes 

N-l 



k=l 



+ u-S = 0. 

(6.25) 

Here /i^ = /i e*^'^'"'/^. Evaluating the sum, one arrives at an algebraic equation 

{l+u + \nfx^)S = /i (^1 + ln(/i5) j . (6.26) 

We can solve Eq. fl6.26p numerically. The result is summarized in Fig. [7] which shows 
the dependence of S and |np on fi. The latter dependence can be inferred from 
Eq. (ESOD, 

— \n\ = In = ln/i5. (6.27) 

Vanishing of |np at certain /i* = m*/A delineates the Higgs and the Cou- 
lomb/confining phases. At that point, a experiences a vertical slope. To the left 
of the phase transition point, |np becomes negative — the analysis of Eq. fl6.26p is 
not valid. Strictly speaking, Higgs phase vacuum equations become invalid as soon 
as I rap reaches zero, one needs to deal with the Coulomb phase. 

Demanding the derivative dS/dfi in Eq. fl6.26p to be infinite, one arrives at an 
equation governing the phase boundary: 

fxl - In/i^ = 1 + u, (6.28) 

together with a useful relation 

fi, = ^. (6.29) 
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Figure 7: Dependence of a and coupling constant |np on m in the Higgs phase for u = 2. 
The former dependence is obtained from Eq. (I6.26p . while the latter is then reconstructed 
from Eq. IKTlh . 



u 




Figure 8: The phase diagram of the twisted-mass deformed heterotic CP(A^ — 1) theory. 
Variable u denotes the amount of deformation, u = ^|a;p. The phase transition boundaries 
are determined by Eq. ()6.28p . 
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The first equation has two solutions, the corresponding curve shown on Fig. [HI The 
larger solution determines the Coulomb / confining-Higgs phase boundary. Its asymp- 
totics at large u, //^ ^ In/i^, is given by (c/. (16.141) ) 



fx^:H — Vl + M . (6.30) 

The smaller solution, as was found in Section 16.4. determines the location of the 
Strong- Coulomb phase boundary. Therefore, single equation (I6.28P defines the whole 
phase diagram of the theory, see Fig. [SI When the deformation parameter u vanishes, 
the whole Coulomb/confining phase shrinks to a point m = A. 

We can now use Eq. (I6.26P to derive a vacuum equation for the coupling constant 
2/3rcn = I'^P- From Eq. (I6.27P one can see that the coupling constant vanishes at 
the point (16. 29 p . confirming once again that a phase transition takes place. Plugging 
Eq. (I6.27P into Eq. (I6.26p . one arrives at 

At large u, one recovers Eq. (16.231) . and therefore the same analysis of the solutions 
sketched after Eq. (I6.23P applies to Eq. (I6.3ip — there are two solutions, only one 
of which is physical. This is the solution shown in Fig. [71 

We now turn to the question of the energy density in the Higgs phase. In terms 
of the dimensionless variables, expression (15. ip can be brought into the form 



Air 



- {S - + uS^. 

Evaluating the sums for large TV, and using vacuum equation (I6.26p . one comes to 
the expression 

= (/i^ - S') In/.^ - Sf - uS\ (6.33) 

It is now straightforward to see that this energy density matches the Coulomb 
phase energy density at the phase transition point /i,,. One eliminates S from 
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Figure 9: Energy densities for Coulomb/confining and Higgs phases for u = 2. Conjunction 
occurs exactly at fj, = fj,^:. 



Eq. (16.331) using relation f l6.29p . and resolves the logarithm In/i^ via (16 .28 p . The 
result is 



iV A2 A2 
Figure [9] illustrates this. 



/i^ - (2 + u)fit + 1 . (6.34) 



7 More on the Coulomb/confining phase 

As was shown above (Sects. 16.11 and 16. 3p . both in the strong coupling and Higgs 
phases the Z2N symmetry is spontaneously broken down to Z2 while in the Cou- 
lomb/confining phase this symmetry remains unbroken (see Sect. 16. 2p . In the former 
two phases we have A^ degenerate vacua, while in the later phase the theory has a 
single vacuum. Just like in non-supersymmetric CP(A^ — 1) model |32] the vacua split, 
and A^ — 1 would-be vacua become quasivacua, see [10] for a review. The vacuum 
splitting can be understood as a manifestation of the Coulomb/confining linear po- 
tential between the kinks [311 [6] that interpolate between the true vacuum, and say, 
the lowest quasivacuum. The force is attractive in the kink-antikink pairs leading to 
formation of weakly coupled bound states (weak coupling is the manifestation of the 
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k=0 k=l k=0 



Figure 10: Linear confinement of the kink-antikink pair. The sohd straight hne repre- 
sents the ground state. The dashed hne shows the vacuum energy density of the lowest 
quasi vacuum. 

1/N suppression of the confining potential, see below). The charged kinks (i.e. the n 
quanta) are eliminated from the spectrum. This is the reason why the n fields were 
called "quarks" by Witten |6]. The spectrum of the theory consists of nn- "mesons." 
The picture of confinement of n's is shown in Fig. [TDl 

This Coulomb / confining linear potential can appear only provided that the pho- 
ton remains massless, which is certainly true in the pure bosonic CP(A^ — 1) model. 
As was pointed out by Witten [6], in the supersymmetric M = (2, 2) CP(A^ — 1) model 
the photon acquires a nonvanishing mass due to a chiral coupling to fermions. In 
particular, with vanishing twisted masses the photon is massive in both = (2, 2) 
[0] and M = (0, 2) [H]. Below we will calculate the photon mass in the model f l3.6p 
and show that it does vanish in the Coulomb/confining phase considered in Sect. 16.21 
This is in accord with the unbroken Zn symmetry detected in this phase. It guar- 
antees self-consistency of the picture. 

To this end we start from the one-loop effective action which is a function of 
fields from the gauge supermultiplet [A^, cr and A). After integration over ra* and C,^ 
in the strong coupling or Coulomb/confining phases (at n = 0) the bosonic part of 
this effective action takes the form [9] 

5eff = J £x!^^F^^ + ^\d^a\^ + V{a) + V2ib5a-b5a)F*Y (7.1) 

where F* is the dual gauge field strength, 

F* = \e,.F,, , (7.2) 

while V{a) can be obtained from (15.11) by eliminating D by virtue of its equation of 
motion (15. Sp . This was done in the closed form for m = in [9], see also Sec. 18.21 
Here and and h are the coupling constants which determine the wave function 
renormalization for the photon, a and sigma-photon mixing respectively {5a is the 
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Figure 11: The wave function renormalization for a. 



quantum fluctuation of the field a around its VEV). These couphngs are given by 
one-loop graphs which we will consider below. In the m = case these graphs were 
calculated in [9]. 

The wave-function renormalizations of the fields from the gauge supermultiplet 
are, in principle, momentum- dependent. We calculate them below in the low-energy 
limit assuming the external momenta to be small. The wave-function renormalization 
for a is given by the graph in Fig. [TT] A straightforward calculation yields 



K=^y^^ — • (7-3) 

The above graph is given by the integral over the momenta of the ^ fermions propa- 
gating in the loop. The integral is saturated at momenta of the order of the ^ mass 
\\pio- + mi|. 

The wave function renormalization for the gauge field was calculated by Witten 
in [6] for zero masses. The generalization to the case of nonzero masses takes the 
form 



1 1 

- = — T 



(7.4) 



.3 + Iv^o- + mj|2 3|v/2o- + mi| 

The right-hand side in Eq. f l7.4p is given by two graphs in Fig. [T21 with bosons it 
and fermions in the loops. The first term in fl7.4p comes from bosons while the 
second one is due to fermions. 

We see that although both Ay, and a were introduced in (13. 6p as auxiliary fields 
(in the limit Cg — )■ o6) after renormalization both couplings, and e^, become finite. 
This makes these fields physical 
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The (Scr) F* mixing was calculated by Witten in [6] in the massless N" = (2,2) 
theory. This mixing is due to the chiral fermion couplings which can make the photon 
massive in two dimensions. In the effective action this term is represented by the 
mixing of the gauge field with the fluctuation of a. It is given by the graph in Fig. [131 
Direct calculation gives for the coupling b 

1 V2a + mi 



a + rrii 



With the given value of b the photon mass is 

mph = e^e^m . (7.6) 

We see that in the strong coupling phase, with a ^ 0, the photon mass does not 
vanish. However, in the Coulomb/confining phase in which the VEV of a vanishes, 
6 = 0, and the photon is massless, as expected, 

mph = 0. (7.7) 

The above circumstance can be readily understood from symmetry considera- 
tions too. The field Sa transforms nontrivially under Z2N symmetry; therefore, the 
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nonvanishing value of b in (I7.ip breaks this symmetry. However, Z2N symmetry is 
restored in the Coulomb/confining phase. Hence, b should be zero in this phase, and 
it is. 

If the photon mass vanishes, generating a linear potential, it is not difficult to 
find the splittings between the quasivacua; they are determined by the value of e^. 
From (17.41) we estimate 

E,^,-E,r^e'^ = -\ (7.8) 

We see that the splitting is a effect. This was expected. 



8 Related issues 

In this section we address a few questions which are not necessarily confined to the 
large- limit. Rather, we focus on some general features of our results, with the 
intention to provide some useful clarifications/illustrations. 



8.1 Remarks on the mirror representation 



br the heterotic 



CP(1) in the Umit of small deformation l| 

In this section we will set all twisted masses to zero. The geometric representation 
of the heterotic TV = (0, 2) CP(1) model is as follows [13]: 

-^^hctcrotic = CliidL Cr + [7(M) Cr R {i dL(p^)ijR + H.C.] (8.1) 



- (0t dR <P) + i^R (0t 9, 0) ] - ^t^ 

where the field appearing in the first line is the spinor field on C, a necessary 
ingredient of the A/" = (0,2) deformation [12]. Here G is the metric, R is the Ricci 



''Subscript (m) of J(m) will indicate that we are working in Minkowski space in this section. 
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tensor and x= 1 + 00^; 

2 2 



cf. Eq. f lKG]) . 



We assume the deformation parameter 7 to be small (it is dimensionless) and 
work to the leading order in 7, neglecting 0(7^) effects in the superpotential. The 
kinetic terms of the CP(1) fields and ip contain in the normalization while 7 in 
the first line is defined in conjunction with the Ricci tensor, so that there is no in 
front of this term. This convention is important for what follows. 

Now, let us remember that the undeformed A/" = (2, 2) CP(1) model has a mirror 
representation [IOl[TT], a Wess-Zumino model with the superpotential 

W,nirror = ^ + ^) ' (S'^) 

where A is the dynamical scale of the CP(1) model. The question is: "what is the 
mirror representation of the deformed model (18. ip . to the leading order in 7?" 

Surprisingly, this question has a very simple answer. To find the answer let us 
observe that the term of the first order in 7 in (18. ip is nothing but the superconformal 
anomaly in the unperturbed A/" = (2, 2) model (it is sufficient to consider this anomaly 
in the unperturbed model since we are after the leading term in 7 in the mirror 
representation). More exactly, in the Af = (2, 2) CP(1) model [l3l 



7,J'a = -^R{du<P^)iYi^)a, (8.4) 

where J'^ is the supercurrent. In what follows, for simplicity, numerical factors like 
2 or TT will be omitted. Equation (18. 4p implies that the 0(7) deformation term in 
(18.11) can be written as 

A£ = 7(„)Ci?(7M'/^)L (8.5) 

Since (18. 4p has a geometric meaning we can readily rewrite this term in the mirror 
representation in terms of Wmirror. Indeed, in the generalized M = (0, 2) Wess- 
Zumino model the term proportional to 7(A/)Cii is ^5] 



AC = CR^pLn', H' = m/dY, (8.6) 
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where "H is the /i-superpotential|f| Moreover, 

(7^^^)l = (WVl)™ + 0(7). (8.7) 
Substituting Eq. (18.70 in (I8.5P and comparing with (18.60 we conclude that 

n = 7(M) . (8.8) 

In principle, one could have added a constant on the right-hand side, but this would 
ruin the Z2 symmetry inherent to the M = (0,2) CP(1) Lagrangian. The constant 
must be set at zero. The scalar potential of the JV = (0, 2) mirror Wess-Zumino 
model is Il5] 



V = |Wf + = + |7,M)r . (8.9) 

where Wmirror is given in (18.31) . The second equality here is valid in the small- 
deformation limit. 

At 7 7^ it is obvious that V > and supersymmetry is broken. The Z2 
symmetry apparent in (18.90 is spontaneously broken too: we have two degenerate 
vacua. 

8.2 Different effective Lagrangians 

In this section we will comment on the relation between the effective Lagrangian 
derived in Sect. [5] from the large- iV expansion and the Veneziano-Yankielowicz effec- 
tive Lagrangian based on anomalies and supersymmetry. For simplicity we will set 
mj = in this section. Generalization to mj 7^ is straightforward. We assume the 
heterotic deformation to be small, m <^ L 

The expansion allows one to derive an honest-to-god effective Lagrangian 
for the field a, valid both in its kinetic and potential parts. The leading order in 
1/N in the potential part is determined by the diagram depicted in Fig. [TT] which 
gives at zero twisted masses 

/:H„ = ^^|9,ap, (8.10) 



^It is worth noting that in [45] the h superpotential H was denoted by S. Note that a broad 
class of the (0, 2) Landau-Ginzburg models were analyzed, from various perspectives, in [^HTIHS] . 
The prime interest of these studies was the flow of the (0,2) Landau-Ginzburg models to non-trivial 
(0,2) supcrconformal field theories [31127]; ^-nd Af — (0,2) analogs of the topological rings in the 
N = (2, 2) theories [3H]- 
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see (I7.3p . The virtual momenta saturating the loop integral are of the order of the 
^ mass \/2|o"|. Up to a numerical coefficient this result is obvious since the field a 
has mass-dimension 1. 

The potential part following from calculations in Sect. [5] is 



2kl 



, 2|a|2 ^ ■ 
In -1+u 



(8.11) 



All corrections to flS.lOp and flS.lip are suppressed by powers of For what 

follows it is convenient to recall the dimensionless variable S (see Eq. (16. lip ). 



S 



a 



A 



Then the large- A^ effective Lagrangian of the a field takes the form 



A^ 



^e.ff — — 



4.{y^'^'^'^'' + ^' [l + |5nin|5P-l + .)]}. 



(8.12) 



(8.13) 



On the other hand, the Veneziano-Yankielowicz method |12] produces an effec- 
tive Lagrangian in the Pickwick sense. It realizes, in a superpotential, the anomalous 
Ward identities of the underlying theory and other symmetries, such as supersymme- 
try, and gives no information on the kinetic part. We hasten to add, though, that in 
two dimensions in the undeformed[^ CP(A^ — 1) models the Veneziano-Yankielowicz 
superpotential Ww = S In S (for twisted superfields) obtained in [35], |36l [7] happens 
to be exact. This was mentioned above more than once. 

In terms of the scalar potential for the a field the Veneziano-Yankielowicz con- 
struction has the form 



yvY = ^ 



N , V2a 
— In — : — 
27r A 



—u2\a\^ 
in 



(8.14) 



The kinetic term (that's where comes from) was not determined; however, we can 
take it in the form obtained in the large- A^ expansion, see (I8.10p . since it is scale 
invariant and, hence, does not violate Ward identities. 
Combining 

el = ^2\a\' (8.15) 



^Thc key word here is "undeformed" , i.e. with no heterotic deformation. 
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(see [9]) with (18 .14^ we arrive at 



In 



cr 



A 



2\aYu 



= ^{j^,\d,S\' + ^'m'\\nSf + \S\'u]Y (8.16) 

It is obvious that the potential in fl8.13p is drastically different from that in (18.161) . 
For instance, (I8.13P contains a single log, while (I8.16P has the square of this logarithm. 
We will comment on the difference and the reasons for its appearance later. Now, let 
us have a closer look at the minima of fl8.13p and fl8.16p . The variable S is complex, 
and there are A^ solutions which differ by the phase. 



= |iS^,| exp 



2TTk 

1^ 



k = 0,1,...,A^- 1, 



(8.17) 



for a more detailed discussion see Ref. Each of these solutions represents one of 
the A^ equivalent vacua. This feature is well-known, and we will omit the phase by 
setting k = 0. Thus, we focus on a real solution. The minimum of (18.131) lies at 



5* = e' 

while the corresponding value of Vcs is 

A^ 



-u/2 



Kff(5*) = — A2(l-e-") 
Arc ' 



.19) 



At the same time, the minimum of (I8.16P lies at 

1 



exp 



2+V4 



u 
2 



■u/2 



U 



+ 



(8.20) 



implying that 



Vvy(5,) = ^ A^ (1 - ^fT^) exp (-1 + ^JY^^) 
47r 



47r ^ 



e-)(l-^ + ...). 



(8.21) 
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The a masses are 




4A2e-"(l-n), 
4A2e-"(l-n)(l-n2 + ...), 



(8.22) 



for ( I8.13P and fl8.16p respectively. The positions of the minima, the a masses as well 
as the vacuum energy densities in these two cases differ by 0{u^) in relative units. 
They coincide in the leading and next-to-leading orders in u, however. 

There are two questions to be discussed: (i) why the effective Lagrangians f l8.13p 
and fl8.16p . being essentially different, predict identical vacuum parameters in the 
leading and next-to-leading order in u; and (ii) why the parameters extracted from 
the 1/N and Veneziano-Yankielowicz Lagrangians diverge from each other at O(u^) 
and higher orders. 

The answer to the first question can be found in [3S]. While the 1/N Lagrangian 
is defined unambiguously, the Veneziano-Yankielowicz method determines only the 
superpotential part of the action. The kinetic part remains ambiguous. We got used 
to the fact that variations of the kinetic part affect only terms with derivatives, which 
are totally irrelevant for the potential part. This is not the case in supersymmetry. 
The correct statement is that variations of the kinetic part term, in addition to 
derivative terms, contains terms with FF, which vanish in the vacuum {F = 0) but 
alter the form of the potential outside the vacuum points (minima of the potential). 
The only requirement to the kinetic term is that it should obey all Ward identities 
(including anomalous) of the underlying microscopic theory. For instance, in the 
case at hand, the simplest choice InE InS does the job. However, 



which reduces to 1 in the points where F = and changes the expression for F (and, 
hence, the scalar potential) outside minima (i.e. at F ^ 0). 

The answer to the second question is even more evident. The Veneziano-Yankie- 
lowicz Lagrangian f l8.16p reflects the Ward identities of the unperturbed CP(A^ — 1) 
model. That's the reason why the predictions following from this Lagrangian fail at 
the level O(m^), but are valid at the level 0{u). We remind the reader that it was 



In S In S 1 + 



'2lnS) (D^lnS) 



does the job as well. In this latter case there is an additional factor 



[l + FF/(aV) + ...] 
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Figure 14: The kinks are represented by the n fields at |mj| < A. 

shown in [13] that the vacuum energy density at the level 0{u) is determined by the 
bifermion condensate in the conventional (unperturbed) CP(A^ — 1) model. 

One last remark is in order here. The kinetic term flS.lOp is not canonic and 
singular at a = 0, implying that this point should be analyzed separately. One 
can readily cast f lS.lOp in the canonic form by a change of variables. Upon this 
transformation cr — t- 5" = (assuming for simplicity a to be real and 

positive), the transformed potential (18. lip develops an extremum at a = (i.e. 
a — )■ — oo). This extremum is maximum rather than minimum. Indeed, at n = 

^pot = ^ (^ - 1) + const. (8.23) 

It is curious to note that (18.23P exactly coincides with the (two-dimensional) dilaton 
effective Lagrangian derived in [51] on the basis of the most general (anomalous) 
scale Ward identities. 

8.3 When the n fields can be considered as solitons 

Long ago Witten showed [6] that the n fields in fact describe kinks interpolating 
between two neighboring vacua picked up from the set of degenerate supersym- 
metric vacua of the A/" = (2,2) sigma model. The above statement refers to the model 
with no twisted masses. (See Fig. [H]). Here we will discuss the physical status of 
these states, and the BPS spectrum at large, as the twisted mass parameter evolves 
towards large values, |m|/A ^ 1. will be assumed to be large so that we can use 
the large- iV solutions. 

In the undeformed A/" = (2,2) theory two distinct regimes are known to exist. 
At large m the theory is in the Higgs regime, while at small m it is in the strong 
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coupling regime. There is no phase transition between the two regimes, since the 
global Zn symmetry of the model is spontaneously broken in both. An apparent 
discontinuity of, say, the derivative of the vacuum expectation of cr at m = A (see 
Fig. [1]) is an artifact of the N oo limit. However, the BPS spectrum experiences 
a drastic change in passing from small to large m. In particular, at |m| < A, the 
masses Mi of the fields in the vacuum |0) are 

= \A - TJiif , i = 0,1,...,N - 1 . (8.24) 

This follows from Eq. (13.61) . The mass degeneracy of the kink iV-plet is gone since 
the twisted mass terms (12.21) break the SU(iV) symmetry of the model leaving intact 
only U(l)^~^ and Zjy-symmetries. 

On the other hand, at large m, in the Higgs regime at weak coupling, the n fields 
no longer describe solitons; rather they represent elementary excitations. In each 
of the vacua there are 2(A^ — 1) real elementary excitations. Say, consider the 
vacuum in which develops a VEV. The phase of n° is eaten by the Higgs mech- 
anism providing the mass to the photon field. The modulus of is excluded from 
dynamics by the constraint \n^f = 2/3. The elementary excitations are described by 
n^, n^, n^~^. The same equation (13. 6p implies that the masses of these elemen- 
tary excitations are0 (at large A^), 

= \mo - rrikf = (m^^ , A; = 1, 2, AT - 1 . (8.25) 

since the vacuum expectation of a in the vacuum under consideration is a ~ mg = m. 
For the values of k that do not scale with A^, each such mass in (18.251) scales as 1/A^ 
at large A^. This is the consequence of our ZAr-symmetric choice of the twisted mass 
parameters (12.21) . 

What about the kink masses in the Higgs regime? (For brevity we will refer to 
them as the weak-coupling regime (WCR) kinks, although this is not quite precise. 
Indeed, if m is in the Higgs domain but ~ a few x A, the coupling constant is of order 
1.) The masses of the WCR kinks can be found from the Veneziano-Yankielowicz 
effective superpotential [37], which is exact in the model under consideration, 

= — > V2(7-mJ In^^ V2a. 8.26 

47r ^ V / A An ^ ' 

1=0 

^*^In fact, equation Mk = |?7io ^ '>^k\ is exact. 
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For definiteness is assumed to be even. One should be careful with the logarithmic 
function. It is obvious that fl8.26p is defined up to imp/ 2 times any integer number. 
For the kinks interpolating between the vacua |0) and \p) one has [37] 



Mkink = 2 |AW| = 2 



WvY(fTo) - >VvY(ap 



- (mo - nip) q 



(8.27) 



where the logarithmic function in fl8.26p is defined with the cut along the negative 
semi-axis, and q is an arbitrary integer. Since the vacuum values of sigma satisfy 



N 



N 



and |m/A| > 1, we have \/2ap = nip, with the exponential accuracy. If ao = mo = m 
and \p\ < N/i, in calculating >Vvy(o"o, (Xp) we do not touch the cut of the logarithm, 
and then we can use the following expression!"] 



Wvy(a) 



47r \ eA 



N 



4ttN 



V2aFj,{y) 



FN{y) 



k=l 



N 



k{kN-l) 



m 



V2 



a 



N 



f8.28) 



Note that is finite at — oo and \y\ < 1. As a result, at \p\ < N/A we arrive at 



M, 



kink 



N ni nlpl 
sm 



71 N 

The lightest mass is obtained by setting q = 



, m 

eA 



N 



kink* 



A^m 7r\p\ 
sm 



, m 

eA 



(8.29) 



(8.30) 



TV N 

At p ~ A^*^ and |m/A| ^ 1 the kink mass does not scale with N. In addition, it is 
enhanced by the large logarithm Inm/eA ~ [Ng"^)'^ compared to the masses of the 
elementary excitations. If p ~ A^^, the scaling law is 

Mkink ~iVrnln(m/A), 



-'^^One might be tempted to use Eq. (|6.6p here, but this does not work since in this expression we 
have a different argument of the logarithm. 
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cf. Eq. (I4.16p . Each kink has a tower of excitations on top of it, corresponding 
to g 7^ 0. The existence of excitations is due to the fact that in addition to the 
topological charge, the kinks can have U(l) charges [28]. Excitations decay on the 
curve of marginal stability (CMS) |m| = eA, see |53j . 
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9 Conclusions 



In this paper we presented the large-A^ solution of the two-dimensional heterotic 
N = (0,2) CP(A^ — 1) model. Our studies were motivated by the fact that this 
model emerges on the world sheet of non-Abelian strings supported in a class of 
four-dimensional M = 1 Yang-Mills theories. The non-trivial dynamics which we 
observed - with three distinct phases, confinement and no confinement, and two 
phase transitions - must somehow reflect dynamics of appropriate four-dimensional 
theories. If so, we open a window to a multitude of unexplored dynamical scenarios 
in A/" = 1 theories. But this is a topic for a separate investigation. 

The heterotic CP models in two dimensions are of great interest on their own. 
For instance, heterotically deformed weighted CP models with twisted masses exhibit 
even a richer phase diagram with highly nontrivial phases. The study of these models 
in the large- limit is in full swing now [5^ . 
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A Notation in Euclidean Space 



Since CP(A^ — 1) sigma model can be obtained as a dimensional reduction from four- 
dimensional theory, we present first our four-dimensional notations. The indices of 
four-dimensional spinors are raised and lowered by the SU(2) metric tensor, 



where 



ea/3 



r n 

-1 



and 





1 



The contractions of the spinor indices are short-handed as 
The sigma matrices for the euchdean space we take as 



I r 



(A.l) 



(A.2) 



(A.3) 



(A.4) 



where are the Pauli matrices. 

Reduction to two dimensions can be conveniently done by picking out x° and as 
the world sheet (or "longitudinal" ) coordinates, and integrating over the orthogonal 
coordinates. The two-dimensional derivatives are then defined to be 



d 



'R 



do + id; 



'3 ! 



Or 



do - ids. 



(A.5) 



One then identifies the lower-index spinors as the two-dimensional left- and right- 
handed chiral spinors 



— 6 ) 



a = 6 



(A.6) 



With these assignments, the dimensional reduction for the contracted spinors then 
takes the following form 



^aA" = -eiijAi], = ^[rXl]. (A.7) 

For two-dimensional variables, the CP(A^ — 1) indices are written as upper ones 
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and as lower ones for the conjugate moduli 

where / = 1, TV. In the geometric formulation of CP(A^ "1)) global indices are 
written upstairs in both cases, only for the conjugate variables the indices with bars 
are used 

<P\ ?, ^r^ = l,...,iV-l, 

and the metric gij is used to contract them. 

B Minkowski versus Euclidean formulation 

In the bulk of the paper we use both, Minkowski and Euclidean conventions. It is 
useful to summarize the transition rules. If the Minkowski coordinates are 

xl, = {t,z], (B.l) 

the passage to the Euclidean space requires 

t^-ir, (B.2) 

and the Euclidean coordinates are 

xl, = {r, z} . (B.3) 

The derivatives are defined as follows: 

= dt + d,, d^ = dt-d,, 
df = dr-zd,, di = dr + zd,. (B.4) 

The Dirac spinor is 

-(::) 

In passing to the Euclidean space = however, \1/ is transformed. 
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Moreover, and are not related by the complex conjugation operation. They 
become independent variables. The fermion gamma matrices are defined as 

< = {1,-^3}, = {1,2^3}. (B.7) 

Finally, 

CE = -CM{t = -iT,...). (B.8) 
With this notation, formally, the fermion kinetic terms in Ce and Cm coincide. 

C Parameters of heterotic deformation 

In this section we list the definitions of parameters of heterotic deformations used in 
previous papers, and their relations to each other, including relations between the 
corresponding parameters in Minkowski and Euclidean spaces. 

It is reasonable to first collect the constants and couplings which necessarily 
accompany the heterotic deformation in theories in which the JV = (0, 2) CP(A^ — 1) 
model arises. If the CP(A^ — 1) model is obtained from a four-dimensional bulk 
theory, then its coupling will be related to the non-abelian gauge coupling of the 
latter theory 

2/3 = -2- (C.l) 

Prom the point of view of the two-dimensional theory, /3 can be understood as the 
radius of the CP(A^ — 1) space 

2 2 / \ 

2/3 = K = -^In l + J^'i'^^M. (C.2) 

^° ^° V ^J=^ J 

The latter Kahler potential induces the round Pubini-Study metric 

G-ii = =^ . (C3) 

Table IC.ll displays the equivalent definitions of the heterotic deformations (the 
bifermionic terms). The subscripts relate the corresponding terms to Minkowski or 
Euclidean space. Equivalence in most cases means equality, except for the relations 
between the Minkowski and Euclidean space expressions, for which the Lagrangians 
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are defined differently, see Appendix |Bl For reasons explained below, different nor- 
malizations of the fields involved in the bifermionic coupling — n\ and C,r — were 
used in definitions of different parameters. In the table, the fields that are normalized 
canonically are shown in bold face, otherwise the fields [i.e. their kinetic terms) are 
normalized to 2/3. 

parameter deformation term 0{N) scaling defined in 



6 




0(1) 


[T2l[T3] 


7(M) 






m 


1(E) 


7(B) = 7(M) 


Vn 




7(B) 


2/3 ■ 7(B) {idLn)iRCR 


0(1) 


m 




2i u XlCr 




m 


u 


N 112 


0(1) 





Table C.l: Parameters of heterotic deformation of the CP(iV — 1) theory. 

One useful normalization (this is not the normalization used in this paper) for 
the supertranslational and superorientational moduli is 2/3, 

£ = 2/3 [|9np + ^RidL^R + CriOlCr + ...] • (C.4) 

In this normalization one has |np = 1. In particular, this naturally comes out when 
the two-dimensional sigma model is obtained from a four-dimensional bulk theory. 
This normalization is useful for studying the sigma model classically. The A/" = (0,2) 
deformation of the sigma model in the gauge formulation is done via the A/" = (0, 2) 
superpotential 

2 

W(a) = 5^. (C.5) 

The latter formula is schematical only, as we do not use N" = (0,2) superfields in 
this paper. Rather, the superpotential leads to the couplings of the type 

2/3 [2z6XlCr + ^iHrXl + ...] . (C.6) 
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The deformation parameter 7 was defined in the |np = 1 normahzation, 

£3 2/3 + (RtdLCR + ... + l,EA^dLn)^RCR] (C.7) 

and is related to the above superpotential deformation 6 as, 

Another form of the bifermionic mixing term is found in CP(1) model, where it can 
be written in terms of real variables S"" and x", see 



^ = P [lid.S^' + Ixr^OlXr + ^CriBlCr + l.E.XRiidLS^'XR + ...] . 

(C.9) 

The relation between S"' and n' (also and depends on the normalization of the 
latter, and in the \n\'^ = 1 case takes the form 

S"" = nr'^n, xl r = nT''^L,R + lL,Rr^n. (C.IO) 



Note the extra factors of 1/2 and 2 required for precise matching between Eqs. (jC.911 
and flOTl) . 

The 2/3 normalization of the kinetic terms also arises in the geometric formulation 
of the heterotic CP(A^ — 1) model. In the latter case, however, one would argue, that 
the supertranslational variable (r has nothing to do with geometry, and therefore 
need not have 2/3 in front of its kinetic term. Parameter 7 was originally introduced 
in [13] in the geometric formulation, in Minkowski space, 

+ 7(M) 9o Cr GijitdL^Wn + .... (C.ll) 

Here, again, the bold face of the variable (^n indicates that its kinetic term is normal- 
ized canonically, whereas the kinetic terms of the orientational variables are naturally 
supplied with a factor of 2/3 coming from the metric Gij. The correspondence between 
the Minkowski and Euclidean space Lagrangians gives 

7{M) = -iliE), (C.12) 
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with 

Gauge formulation appears to be the most useful for studying the quantum effects 
of the sigma model. In the quantum theory, one prefers to depart from |np being 
equal to unity, as the former determines the coupling of the theory and may, in 
principle, vanish. One absorbs all outstanding factors of 2/3 into the definition of n\ 
and C,R. This is the normalization used in the bulk of this paper, 

C = |VXl' + ?z^^^V^^' + IriOUr + i^dLri)iRCR + ••• , (C.14) 



where 7(b) is still related to 5 via Eq. flC.81) . The N = (0, 2) superpotential couplings 
( 1C.6I) become 

2iioXLCR + 2iuJCRXL + (C.15) 

where 

u = V^^- (C.16) 
In the large- limit, a more appropriate parameter appears to be u, 

u = ^|c.r, (C.17) 

which enters the effective potential as 

Kff D ^■u\a\\ (C.18) 

Note, that relation ( 10.161) is classical, while quantum mechanically, all couplings 
run. In the one-loop approximation of this paper, the couplings u and /3 diagonalize 
the renormalization group equations, and are, therefore, the appropriate parameters 
in the quantum regime. In the large- limit, again, parameter u is more preferrable. 

Finally, in Ref. [I3], parameter a is used, which relates to 5 and 7 via 



The relation between the heterotic parameters in Euclidean and Minkowski spaces 
IS given byEq. flOl2ll . Every wh ere in the paper where there is no menace of confusion 
we omit the super/subscripts M,E. 
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D Global symmetries of the CP (AT — 1) model 
with ^iv-symmetric twisted masses0 

In the absence of the twisted masses the model is SU(A^) symmetric. The twisted 
masses (12.21) exphcitly break this symmetry of the Lagrangian (13. 6p down to U(l)^~^, 

n' ^ e^"^n^ eWe^^^e^ el^e*"^el, ^=l,2,...,iV, 

(J ^ (7 , \r,l^\r,l- (D.l) 

where are constant phases different for different i. 

Next, there is a global vectorial U(l) symmetry which rotates all fermions in 
one and the same way, leaving the boson fields intact, 

iR ^ e^^^R. iL^^'^iL. ^ = l,2,...,iV, 

\r ^ e-'^\R. Xl^c-'^Xl, 

^ n\ a . (D.2) 

Finally, there is a discrete Z2N symmetry which is of most importance for our 
purposes. Indeed, let us start from the axial U(l)ij transformation which would be a 
symmetry of the classical action at m = (it is anomalous, though, under quantum 
corrections) , 

iR ^ e^'iR. ei^e-^^ei, ^=1,2,...,7V, 

Xr ^ c'^Xr, Xl^c-'-'Xl, a^e^^V, 

^ n'. (D.3) 

With m switched on and the chiral anomaly included, this transformation is no 
longer the symmetry of the model. However, a discrete Z2N subgroup survives both 
the inclusion of anomaly and m 7^ 0. This subgroup corresponds to 

7fc = ^, fc = l,2,...,iV. (D.4) 

with the simultaneous shift 

e^e-k. (D.5) 
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Sec also the Appendix in Ref. [49] 
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In other words, 



n^-\ (D.6) 

This Z2N symmetry rehes on the particular choice of masses given in (12. 2p . 

When we switch on the heterotic deformation, the transformations ( 1D.6P 
must be supplemented by 

U^e-'^^CR- (D.7) 

The symmetry of the Lagrangian fl3.14p remains intact. 

The order parameters for the symmetry are as follows: (i) the set of the 
vacuum expectation values {{'nP)-, (n^), ■■■ {n'^~^)} and (i) the bifermion condensate 
{^Rj^i}- S^Y' ^ nonvanishing value of (n°) or {^r^ g^f) implies that the Z2N symmetry 
of the action is broken down to Z2. The first order parameter is more convenient for 
detection at large m while the second at small m. 

It is instructive to illustrate the above conclusions in the geometrical formulation 
of the sigma model. Namely, in components the Lagrangian of the model is (for 
simplicity we will consider CP(1); generalization to CP(A^ — 1) is straightforward) 

C CPU) = G[d,^d^<P-\m\'^<f)+'-{^[dnijL + iljl Ol ^Ar) 

-i — (mi)li)R + Thtpl^i)L) 

X 

-^V'ULV^Uij}, (D.8) 



where 



X = l + ^(f), G = ^^. (D.9) 
9o X 



The Z2 transformation corresponding to ( ]D.6|) is 

ij^iPL-^-^li^L. (D.io) 
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The order parameter which can detect breaking/ nonbreaking of the above symmetry 



IS 

m 

7, 



Under the transformation f ID.lOp this order parameter changes sign. In fact, this is 



the central charge of the M = 2 sigma model, including the anomaly |l3l E]. 

E Geometric formulation of the model 

Here we will briefly review the M = (2, 2) supersymmetric CP(A^ — 1) models in the 
geometric formulation. 

E.l Geometric formulation, 7 = 

As usual, we start from the undeformed case. The target space is the — 1- 
dimensional Kahler manifold parametrized by the fields 0*, (^^-^^ i, j = 1, . . . , N — 1, 
which are the lowest components of the chiral and antichiral superfields 

<i>'(x^ + z^7^e), ^^x^'-ieye), (E.l) 

wherein 

xf' = {t,z}, e = e^-i\ v = ^/'V 

7° = 7* = aa , 7^ = 7^ = iai , 75 = 7%^ = t^a • (E.2) 
With no twisted mass the Lagrangian is [IS] (see also 

(E.3) 

where 

G, = (E.4) 
is the Kahler metric, and Rijki is the Riemann tensor [2^ , 

Rfjkfh ~ {GijGkm. + GimGkj) ■ (E-5) 
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In the Euclidean space V' becomes an independent variable. 
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Moreover, 

is the covariant derivative. The Ricci tensor R^j is proportional to the metric 

9^0 



R^ = fNG^. (E.6) 



2 

For the massless CP(A^— 1) model a particular choice of the Kahler potential 

In ( 1 + 

9o 



Km=o = ^ In I 1 + ^ 1 (E.7) 



corresponds to the round Fubini-Study metric. 

Let us briefly remind how one can introduce the twisted mass parameters 
The theory (IE. 30 can be interpreted as an = 1 theory of — 1 chiral superfields in 
four dimensions. The theory possesses — 1 distinct U(l) isometries parametrized 
by t", a = 1, . . . , — 1. The Killing vectors of the isometries can be expressed via 
derivatives of the Killing potentials D"-{(j),(f)), 

^ = -^G^^^ ^ = ^G^^ (E 8) 

dta d(j)^ ' dta d(f)' ' ^ ■ ^ 

This defines the U(l) Killing potentials, up to additive constants. 

The N — 1 isometries are evident from the expression (IE.7P for the Kahler poten- 
tial, 

5<j)' = -t5ta{T%{(t>f , = tSUT-)j^\ a = 1, . . . , iV - 1 , (E.9) 

(together with the similar variation of fermionic fields), where the generators T" have 
a simple diagonal form, 

(ni = «, a = l,...,iV-l. (E.IO) 

The explicit form of the Killing potentials in CP(A^— 1) with the Fubini-Study 
metric is 

D^ = l-^, a = l,...,iV-l. (E.ll) 
61 



Here we use the matrix notation implying that </> is a column 0* and is a row (jy' . 

The isometries allow us to introduce an interaction with — 1 external U(l) 
gauge superfields Va by modifying, in a gauge invariant way, the Kahler potential 

dEZD, 

K„,=o{^,^)^Km{<^,^,V). (E.12) 
For CP(A^ — 1) this modification takes the form 



i^^ = 4ln(l + <le^"^"$) . 



(E.13) 



In every gauge multiplet Va let us retain only the and Ay components of the gauge 
potentials taking them to be just constants, 

Va = -maOil + 75)^ - maOil - 75)^ , (E.14) 

where we introduced complex masses as linear combinations of constant U(l) 
gauge potentials, 

A n . • A n — * 

■m„ = m„ 



rria — Ay + iA'^ , 



1,2, 



(E.15) 

The introduction of the twisted masses does not break A/" = (2, 2) supersymmetry in 
two dimensions. To see this one can note that the mass parameters can be viewed 
as the lowest components of the twisted chiral superfields 

Now we can go back to two dimensions implying that there is no dependence 
on X and y in the chiral fields. It gives us the Lagrangian with the twisted masses 
included |23l|2H]: 



d^e KU^, $, V) = gMN V^'P 1)^0* + t 



(E.I6) 



where G 



didjKm\g=g=o is Kahler metric and summation over M includes. 



besides M - 


= 


= 0,1, 


also M = 


are 










/ 1 





\ 







-1 





gMN = 












u 





-1 0/ 



7^ 



-z(l+75), 7- = z(l-75). (E.17) 
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The gamma-matrices satisfy the following algebra: 



2^? 



MN 



(E.18) 



where the set F differs from F by interchanging of the +, — components, F = F^. 
The gauge covariant derivatives P^'^ are defined as 



(E.19) 



and similarly for T) ip, while the general covariant derivatives D ip are 

D^^' = V^^' + Til ^^0*^ ^' • 



(E.20) 



E.2 Geometric formulation, 7 7^ 

The parameter of the heterotic deformation in the geometric formulation will be 
denoted by 7(^7) (the tilde appears here for historical reasons; perhaps, in the future 
it will be reasonable to omit it; subscript (m) will stress that this section works with 
Minkowski notations). 

To obtain the Lagrangian of the heterotically deformed model we act as follows 
|16j : we start from flE.3l) . add the right-handed spinor field (r, with the same kinetic 



term as in Sect. 13.21 and add the bifermion terms 

7(M) CnGfj ) + 7(M) CBpfj 



9o 
V2 



(E.21) 



Next, we change the four-fermion terms exactly in the same way this was done in 
[13], namely 



1 



R 



ijkl 



ij^^' ) ( v^V'' ) ( ^^V'* ) ( 



L ' 



(E.22) 
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where the first hne represents the last term in Eq. ( 1E.3I) . and we used the iden- 
tity ( 1E.5I) . If one of the twisted masses from the set {mi, m2, ...,171^} vanishes (say, 
0), then this is the end of the story. The masses rria in Eqs. (1E.14I) and (IE. 1511 



m 



N 



are {mi, m2, ...,mAr_i}. 

However, with more general twisted mass sets, for instance, for the Zyv-symmetric 
masses (12.21) . one arrives at a more contrived situation since one should take into 
account an extra contribution. Occurrence of this contribution can be seen [16] in a 
relatively concise manner using the superfield formalism of |13j . 

A>C ~ M j BSl dOn + H.c. , (E.23) 

where is a (dimensionless) M = (0,2) superfield0 

B = {(r (x'^ + le-f^^e) + V29rT^ 9l . (E.24) 

The parameter M appearing in (]E.23P has dimension of mass; in fact, it is propor- 
tional to . 

As a result, the heterotically deformed CP(A^ — 1) Lagrangian with all twisted 
mass parameters included can be written in the following general form: 

C = C( + Cm=o + Cm , (E.25) 

where the notation is self-explanatory. The expression for Cm. is quite cumbersome. 
We will not reproduce it here, referring the interested reader to For convenience, 
we present here the first two terms. 



■^C + Cm=o — Cr Cr 



7(„)-^CH,G.j(z9z.0t^>lj + H.c. 



+f (1 - l7(M,r) {G^l^'Ji^l) {g^^A^^P'r) , (E.26) 



14 



This means that B is the superfield only with respect to the right-handed transformations. 
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where we used ( 1E.5I) . The above Lagrangian is A/" = (0, 2) -supersymmetric at the 
classical level. Supersymmetry is spontaneously broken by nonperturbative effects 
[T2I |9]. Inclusion of Cm spontanesously breaks supersymmetry at the classical level 
(see Eq. f l3J[5|) and Eq. (2.11) in [l6]). 

The relation between 7 and 5 is as follows [16] : 

i 7(Af) = 7(E) = \/2 —====, (E.27) 
a/1 + 2|(^|^ 

implying that 7 does not scale with in the 't Hooft limit. See Appendix [B] for 
details on relation between Euclidean and Minkowski notations. 
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